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Nonlinearity and entanglement arc two important indications by which physical systems can be 
identified as non-classical. We study the dynamics of the resonant interaction of up to N=3 two-level 
systems and a single mode of the electromagnetic field. We focus on a single excitation dynamically 
shared in this fourpartite system. We observe coherent vacuum Rabi oscillations and their nonlinear 
\/N- speed up in frequency by tracking the populations of all three qubits and the resonator. We 
use full quantum state tomography to verify that the dynamics generates a maximally entangled 
3-qubit W-state when the cavity returns to the vacuum state. 



EflFects related to the interaction of light and matter 
can be investigated in many different physical systems. 
The description of this interaction can be simplified to 
the interplay of identical two-level systems and a single 
mode of an electromagnetic field. This theoretical ab- 
straction is known as the Dicke- or Jaynes-Cummings 
model [1, 2]. In recent years it became feasible to inves- 
tigate the complex interaction of light with multiple two- 
level systems as experiments approached a higher level of 
control over collections of individual quantum systems. 
The coupling strength (transition amplitude) of N two- 
level systems and a single mode increases as \fN . Clear 
evidence of this nonlinear behavior was observed in spec- 
troscopic measurement with few atoms [3 6], large en- 
sembles of atoms using cold gases [7-9], ion Coulomb 
crystals [10] and by using superconducting qubits cou- 
pled to a transmission line resonator [11, 12]. 

More insight on the dynamics of collective systems can 
be gained by time resolved measurements of energy ex- 
change between its individual components. When multi- 
ple two-level systems are resonantly coupled to a single 
mode, this process is called collective vacuum Rabi os- 
cillations. The collective coupling strength g defines the 
frequency of these oscillations. Here, we restrict our in- 
vestigation to the initial state in which the cavity is popu- 
lated with exactly one photon. In this case the oscillation 
involves a single photon that is continuously absorbed 
and recmitted by all N two-level systems. In fact each 
two-level system absorbs the photon with equal proba- 
bility 1 /N. Our lack of knowledge about which two- level 
system actually absorbs the photon leads to entangled 
states which are known as W-states [13]. The nonlinear 
enhancement of the coupling strength speeds up the gen- 
eration of W-states when N is increased. The preparation 
of Dicke- and W-states has been demonstrated in NMR- 
[14], photonic- [15], ion trap- [16] and circuit QED [17] 
systems. 

Circuit QED, in which a controllable number of super- 
conducting qubits arc strongly coupled to a single mode 
of a radiation field defined by a transmission line res- 



onator, offers a suitable set of tools to investigate collec- 
tive interactions. So far the resonant collective dynamics 
of a circuit QED system has been observed in time re- 
solved measurements of two phase qubits coupled to a 
transmission line resonator [12]. Time resolved measure- 
ment, however, have not been carried out with more than 
two qubits and clear evidence of entanglement, for ex- 
ample by measuring the density matrix of the generated 
states, is still pending. Generally, the generation and 
observation of multi-qubit entangled states is quite chal- 
lenging in these systems and was achieved only recently 
for up to three qubits [17-19]. 

Here, we report time domain studies of the collective 
oscillation of three fully controllable transmon-type su- 
perconducting qubits [20] coupled to a single electric field 
mode of a microwave resonator. The coupling of multi- 
ple two-level systems to a single mode is described by the 
Tavis-Cummings Hamiltonian [21] 

Htc = huj^a^a+^ (^'^^^^ ^Qjio^^J + ^/a) j • (1) 

u)^ is the frequency of the field mode, uij is the transition 
frequency between ground state \g) and excited state |e) 
of the two-level system j, a and a} are the creation and 
annihilation operators of the mode and are the cor- 
responding operators acting on the qubits. In our cir- 
cuit QED implementation of the Hamiltonian "Htc we 
use the first harmonic mode of a coplanar transmission 
line resonator at Wj./27r ~ 7.023 GHz with quality fac- 
tor Q ~ 14800. An optical microscope image of the 
sample is shown in Fig. 1(a). The resonator is used for 
the resonant exchange of a single excitation and for dis- 
persive joint thrce-qubit readout by measuring its trans- 
mission properties [22]. All qubits are located at the 
antinodes of the first harmonic mode of the resonator. 
The frequency of each qubit is approximately given by 
hiOj{^) « ^y8EcEj{(f)) — Ec where the Josephson en- 
ergy depends periodically on the flux penetrating the 
SQUID loop according to Ej{(t>) = Ejmax\co^{T''4>/(f'o)\- 
The maximum Josephson energies Ej{0)/h for the three 
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FIG. 1. (a) Optical microscope image of the sample with three 
qubits (A,B,C) capacitively coupled to a coplanar waveguide 
resonator, (b) Enlarged view of the transmon qubit C. Each 
qubit is equipped with individual local charge and magnetic 
flux-bias lines. 



qubits are (26.8, 28.1, 25.7) GHz and their charging en- 
ergies Ec/h are (459, 359, 358) MHz. The anharmonic- 
ity of the transmon energy levels depends on Eq and 
is chosen such that the validity of the two-level approx- 
imation is ensured while keeping the charge dispersion 
low. The maximum transition frequency of the qubits 
a;j(0)/27r = (9.58,8.65,8.23) GHz is designed such that 
all three qubits can be tuned into resonance with the res- 
onator. In the steady state the qubits are flux biased at 
Wj/27r ~ (6.11, 4.97, 7.82) GHz using the quasistatic mag- 
netic field generated by three superconducting miniature 
coils positioned underneath the chip, such that |Aj | ^ gj 
with detuning Aj — ujj — . We measure qubit dephas- 
ing times T2 — (100, 140, 440) ns and qubit relaxation 
times Ti = (2.1, 1.8, 1.0) fj,s. Tuning of the transition 
frequencies on the nanosecond timescale is achieved by 
injecting current pulses into on-chip flux control lines (see 
Fig. 1(b)). Both for the coils as well as for the flux gate 
lines we determined the full coupling matrix to compen- 
sate for cross coupling. The coupling strengths gj of the 
qubits (A,B,C) extracted from spectroscopic measure- 
ments are {gA,9B,9c)l'^ = (-105.4,110.8,111.6) MHz. 
The negative coupling constant of qubit A originates from 
the TT-phase difference of the first harmonic mode between 
the center of the resonator and its coupling ports. 

By applying phase controlled truncated Gaussian 



DRAG pulses [23] through the on-chip charge bias lines 
(see Fig. 1(b)) each qubit can be prepared in an arbi- 
trary superposition state tp = a\g) -\- I3\e) and read out 
jointly [22]. The extraction of the N qubit populations 
needs 2^ single qubit rotations and the tomographic ex- 
traction of the full density matrix requires 2^^ pulses. 
The pulse sequence we have implemented to observe the 
one qubit vacuum Rabi oscillation is depicted in Fig. 2 (a, 
left column). First, a 7r-pulse is applied to a single qubit 
far detuned from the resonator, such that the system can 
be described by a product state |e) ig) |0). Then the qubit 
is tuned into resonance with the resonator using a flux 
pulse of variable duration r. To reduce large overshoots 
during resonant interaction the qubit is tuned to an inter- 
mediate buffer frequency before the final resonant pulse 
is applied. On resonance the energy eigenstates of this 
system containing n = 1 excitations and N qubits are 
|n = l,iV = l±) = l/\/2(|g, 1) ± |e,0)) and the initial 
excited qubit state undergoes vacuum Rabi oscillations 
between [g, 1) and |e, 0). The frequency of the oscillations 
is given by = 2a/P~+A2 with VIq = 2g for A = 0. The 
amplitude of the vacuum Rabi oscillation is maximal on 
resonance where the excitation is fully exchanged. After 
the resonant flux pulse of length r the qubit is detuned 
from the resonator again and the energy exchange pro- 
cess is stopped. In addition to the dispersive readout for 
the qubit population we determine the resonator popula- 
tion by measuring the average photon number, as given 
by the time integrated power (a^a) at the output of the 
cavity once the oscillation is stopped [24] . As can be seen 
in Fig. 2 (a, right column) the qubit population is oscil- 
lating with frequency of 112.0 MHz and is exactly out of 
phase with the photon field oscillating at an experimen- 
tally extracted frequency of 111.2 MHz. Both values are 
in good agreement with the spectroscopically obtained 
value. 

We extended the procedure described above to two and 
three qubits as shown in the pulse diagrams in Fig. 2(b) 
and 2(c). We study the collective vacuum Rabi oscilla- 
tions with initially a single excitation in the resonator. 
This initial state is prepared by transferring the full ex- 
citation of the first qubit to the resonator by adjusting 
the interaction time between the qubit and the resonator 
to To = 7r/25. Then the second and third qubit are 
tuned into resonance and the resonant collective interac- 
tion proceeds for time r. All qubit populations oscillate 
simultaneously out of phase with the cavity. To obtain 
all qubit populations by a tomographic state reconstruc- 
tion with a joint readout a set of single qubit rotation 
pulses is applied at the qubit steady state frequencies for 
each interaction time t. When the number N of qubits 
taking part in the resonant interaction is increased, the 
frequency of the oscillation scales with ^/N and the am- 
plitude of the individual qubit population decreases to 
1/N such that normalization is fulfilled. The tripartite 
states between which the oscillations occur are |g, 1) 
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FIG. 2. Dynamics of the collective vacuum Rabi oscillations, (a) Pulse sequence for the single qubit vacuum Rabi oscillation 
(QP=Qubit Preparation, SRI=Single Resonant Interaction, ST=State Tomography). The population P is shown for the 
qubit C (blue) and the cavity (violet), (b) Pulse sequence for the two qubit collective vacuum Rabi oscillation (CRI=Collective 
Resonant Interaction). The population is shown for two qubits A (green) and C (blue), for the ground state (gray) and the 
cavity (violet), (c) Pulse sequence for the three qubit collective vacuum Rabi oscillation. The population is shown for qubit A 
(green), B (brown) and C (blue), for the ground state (gray) and the cavity (violet). 



and l/y2(|g, e,0) - |e,5, 0)) and equivalently \g,g,g,l) 
and l/V3{\g,g,e,Q) + \g,e,g,0)-\e,g,g,0)) for the four- 
partite state, where we have denoted the states in a bi- 
nary order as |C, B, A, Cavity). The opposite phase of 
the \e,g,g,0) state originates from the fact that qubit A 
has a negative coupling constant. The collective oscilla- 
tion transfers energy back and forth between the qubits 
and the resonator, such that the population in the cavity 
is correlated with the population of the qubit \ggg) state. 
Hence it is instructive to measure the photonic part of 
the state as presented before for the single qubit vacuum 
Rabi oscillations (see Fig. 2(b/c, third column)). 

The qubit populations measured after the collective 
vacuum oscillations are shown in Fig. 2(b) and Fig. 2(c, 
second column). The oscillation frequencies for the res- 
onator population for N=2 and N=3 are 161.8 MHz 
= 1.05 • V2 g/2TT and 195.2 MHz = 1.03 • ^3 5/27r, with 
g being the root mean square [25] of the three spectro- 
scopically obtained coupling constants. These oscillation 
frequencies clearly demonstrate the \/]V-nonlinearity of 
the Tavis-Cummings Hamiltonian in time resolved mea- 
surements. While the collective qubit oscillations are 
in good agreement with the expected dynamics, we ob- 



serve a decrease in visibility with increasing N for both 
the qubit and resonator populations. We attribute this 
observation to an imperfect control of the flux pulse 
amplitudes. Although we extracted a flux pulse cross 
coupling matrix, its accuracy is not sufficient to fully 
compensate the flux cross talk and to exactly tune all 
qubits into resonance with the resonator. During the 
collective oscillations the W-state is created at time 
T ~ T:/{2g\/N), i.e. when the cavity state factorizes 
|0) «) 1/V3{\g,g,e) + \g,e,g) - \e,g,g)). We measured 
the three qubit density matrix at this time by applying 
full quantum state tomography [22] . The result is shown 
in Fig. 3(a) and was obtained by applying the pulse se- 
quence shown in figure 2(c). All three qubits were tuned 
out of resonance after only r = 2.9 ns of collective in- 
teraction. If needed the dynamic phase that is picked 
up during the time when the qubits are detuned can be 
corrected by applying single qubit phase gates or by ad- 
justing the phases of the tomography pulses. Here the 
coherent entries of the density matrix have been numeri- 
cally rotated to correspond to the phases in the expected 
state. This allowed us to characterize the measured state 
by defining the fidelity as {pt \ Pppmlpt) with pppm be- 
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FIG. 3. Measured density matrix of the W-state. (a) Obtained for the collective approach as shown in figure 2. (b) Generated 
by the sequential approach. 



ing the rotated density matrix. After applying a max- 
imum likelihood estimation [26] the resulting fidelity is 
78% with respect to the density matrix pt = \'^w) {'^w\ 
and = 1/V^(lff,5,e) + \g,e,g) - \e,g,g)) being the 
ideal state. As already discussed in the context of time 
resolved measurements the dominant limitation to our 
control over the oscillation is flux cross talk. 

For comparison we have also prepared the W-state 
using a sequential method. Also in this approach we 
distribute a single excitation equally between the three 
qubits using resonant interaction. First, we excite 
qubit C and tune it into resonance with the resonator 
for time ri = arcsin ( )/gc- During this interac- 
tion, only 2/3 of the energy is transferred to the res- 
onator. Next we tune qubit B into resonance for a time 
T2 = arcsin ( \/l/2 )/gB, transferring half of the resonator 
energy to qubit B. Finally we let qubit A pick up the re- 
maining third of the energy from the resonator by bring- 
ing it into resonance for a time T3 = arcsin (1) /gA- Using 
a joint readout and tomography we obtain a density ma- 
trix of the W-state with fidelity 91%. 

As expected the preparation time for the consecutive 
attempt is indeed longer than for the collective approach, 
because the latter method exploits the V iV-enhancement. 
The fact that the higher fidelity was obtained by a proce- 
dure where no flux pulses are injected simultaneously af- 
firms that the main limitation of the collective generation 
is due to residual flux crosstalk rather than decoherence. 

In summary, we have experimentally studied collec- 
tive vacuum Rabi oscillations of three superconducting 
qubits coupled to a microwave resonator. We have fully 
resolved the dynamics of the population in the time do- 
main. The v'iV-scaling of the oscillation frequency of 
the collective vacuum Rabi oscillations and the genera- 
tion of multi-qubit entanglement were studied. Moreover 
we compared the collective method to an approach that 
achieves entanglement sequentially. The techniques used 
and developed allowed for the complete examination of 
collective effects by taking advantage of individual con- 
trol over all involved system components. 
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